Polarization Diffusion from Spacetime Uncertainty 
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A model of Lorentz invariant random fluctuations in photon polarization is presented. The effects 
are frequency dependent and affect the polarization of photons as they propagate through space. 
We test for this effect by confronting the model with the latest measurements of polarization of 
Cosmic Microwave Background (CMB) photons. 
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All approaches to the problem of quantum gravity pre- 
dict that the spacetime itself will suffer from quantum un- 
certainty at Planckian scales. This basic idea is the inspi- 
ration behind many attempts to formulate phenomeno- 
logical models of quantum gravitational effects with po- 
tentially observable consequences. To date, much of the 
effort spent on modeling the effect of spacetime fluctua- 
tions has produced Lorentz symmetry violating models. 
However, with constraints on violations of Lorentz sym- 
metry becoming tighter all the time (see e.g. |T) it is 
more important than ever to discover quantum gravity 
phenomenology that respects Lorentz symmetry. 

That such models can exist has been demonstrated in 
[5] and [3] . That work was motivated by the causal set ap- 
proach to quantum gravity but the scheme is quite 
general and does not depend on any details of the under- 
lying theory except that it should be Lorentz invariant. 
The basic idea is that certain dynamical quantities such 
as particle trajectories are subject to minute, quasi-local, 
random fluctuations due to the uncertainty in spacetime 
structure at the Planck scale. The model-building strat- 
egy is straightforward: identify a space of states for the 
system, work out how Lorentz transformations act and 
hence deduce the most general Lorentz invariant diffusion 
process on that space. 

In the case of a massive point particle, the outcome 
of this strategy is an Ornstein-Uhlenbeck-type process in 
which the momentum of the particle undergoes Brownian 
motion on the mass shell in proper time [2] . For massless 
particles, the Lorentz symmetry restricts the process to 
be a one dimensional diffusion in energy - the particles 
always travel on the light cone - but a second indepen- 
dent parameter enters which governs a drift in energy 

In this paper we will apply the strategy described 
above to polarization degrees of freedom as suggested in 
[7] . We model a photon classically as a point particle with 
a spacetime position x^, null momentum /c^ — {k^,k) 
and a polarization state to be identified. A more realis- 
tic description would use wave packets and an even better 
model would take account of the quantal nature of pho- 
tons. For now we assume that this classical state is a 
good approximate description of each of the free stream- 
ing photons produced by astrophysical and cosmological 
sources which reach our detectors. 

The state space for a classical photon is therefore M** x 
T-Lq X B where is 4 dimensional Minkowski spacetime. 



Hq is the 3 dimensional "cone" of future-pointing null 
4- vectors, and B is the space of polarization states which 
we will see is the Bloch sphere. 

The polarization state of a massless particle of momen- 
tum fc'^ can be given by a complex 4-vector a'^ such that 



kt'au 



and a^*ai, 



I. The vector a 



for any complex number A, will describe the same state. 
To eliminate this gauge freedom, we can consider the po- 
larization state to be given by the complex two form, 
P — kAa, whose components, P^,^ — k^a^, — a^k^ satisfy 
the Lorentz invariant conditions 



P^"'k^ =0, 



pt^v*p^ 
pt^u*P 



= 0, 
= 0, 



If fc'' = s'' where s^' := (1,0,0,1) 
components 
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P has the following 
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where Oi and a2 are complex numbers such that |aip + 
|a2p = 1. This corresponds to a polarization vector = 
(0, fli, a2, 0). 

The phase of the 2-d complex unit vector (ai,a2) is 
not relevant for the polarization state of a single photon 
and so the polarization state space has two real dimen- 
sions: it is the Bloch sphere, B = CP^. Let a and /3 be, 
respectively, the usual polar and azimuthal angles on B, 
then they are related to the components of by 



02 =iT^{coa1 



- e'^^ sin - 



e"^^ sin - 



(6) 
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where 7 is an irrelevant phase. The north and south 
poles, a — 0,TT, are the circularly polarized states and 
the equator, a = 7r/2, consists of the linearly polarized 
states. 

Now consider a general photon state (fc'^ , P^^). For a 
general k^, the polarization 2-form Pp,y must be trans- 
formed by a Lorentz transformation that takes k^^ to s^^ 
in order for it to be compared to the standard polariza- 
tion basis and its coordinates on B determined. This can 
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be done using a standard Lorentz transformation defined 
for example in [8]. If P{k) is the polarization 2-form thus 
transformed, then it will have components of the form ([s]) 
and 



P{k), 



(8) 



where = (0, P{k)if), P{k)2o, 0). The {a, fi) coordinates 
of the polarization state on B are then obtained from ^ 
with oi = P{k)io and 02 = P{k)2o . 

In this way, every photon state is specified by coordi- 
nates (a;'', kt", a, /3) on x Hj] x B. 

Under a Lorentz transformation the photon state 
(k^^P^v) transforms in the usual way as a vector and 
2-tensor and it can be shown that this translates into a 
polar rotation on B, a rotation around the north-south 
polar axis generated by Details of these derivations 
will appear elsewhere. 

The Stokes parameters (see e.g. fSl) are a convenient 
way to parameterize the polarization of a beam of electro- 
magnetic radiation. A monochromatic beam with Stokes 
parameters (/, Q, C/, V) can be modeled as a bunch of 
photons with the same momentum k^ and polarization 
states distributed over B. I is the intensity of the beam 
and since our process preserves particle number / is fixed. 
If a beam consists of photons of momentum k^ which are 
all in the same polarization state (a,/?) G B then the 
Stokes parameters of this perfectly polarized beam are 
Q = I sin a cos l3, U = I sin a sin /3 and V — I cos a. 

If the photons have a distribution of polarizations the 
Stokes parameters are weighted by the probability den- 
sity p(a,P) on B e.g. Q — I J^^ sin a cos /3 p{a, (f)) dad/S 
and similarly for U and V. There are many distributions 
that will model a given set of Stokes parameters. For ex- 
ample, an unpolarized beam, Q = U = V = 0, could be 
modeled by a uniform distribution of linearly polarized 
states, or a uniform distribution on the two circularly po- 
larized states alone. In general, the more spread out the 
distribution on B, the smaller the polarization fraction, 

V := v'Q2 + [72 + vyi. 

Having identified the state space of the photon as 
X Hq X B we can deduce the most general Lorentz in- 
variant diffusion process on this space. As described in [3, 
the trajectory in spacetime is simple: the photon moves 
along null lines according to = A:^ where A is affine 
time. Moreover the process on the momentum space Hq 
must not disturb the blackbody nature of the CMBR 
spectrum over the age of the universe. This means that 
we can neglect this effect for the purposes of this paper: 
we assume that the photon's frequency is constant along 
its worldline. We are left with the task of deducing the 
Lorentz invariant diffusion equation on B. 

We refer to coordinates on B as = {a, (3). Then, 
following [TP] , the most general diffusion equation on B 
is 



dp 
dX 



dA{K^''nds{^)~u^p) 



(9) 



where A is affine time, K"^^ is a symmetric, positive semi- 
definite 2-tensor, u'^ is a vector and n is a scalar density 
("density of states") on B. These geometric quantities 
are the phenomenological parameters of the model and 
must be Lorentz invariant. 

There is an embarrassment of choice of parameters be- 
cause Lorentz transformations act as polar rotations only. 
Any tensor, vector or scalar density that does not depend 
on the azimuthal angle, /?, is Lorentz invariant. Free pa- 
rameters that are whole functions do not make for pow- 
erful phenomenology. If, however, we restrict attention 
to the linear polarization states alone (corresponding to 
setting Stokes parameter V to zero), the model recovers 
its predictive power. 

The space of linearly polarized states is the unit circle, 
the equator of the Bloch sphere, and the Lorentz trans- 
formations act as rotations of the circle. The coordinate 
around the circle is f3 and there is, up to a constant factor, 
one Lorentz invariant vector, d /dj3. A Lorentz invariant 
density n must be constant on the circle. We deduce a 
simple diffusion-cum-drift on the circle for the distribu- 
tion p = p{a, j3): 



'a/3: 



(10) 



where c > and d are constants. 

Transforming from affine time to "cosmic time", i.e. 
time in the observatory frame, t — hv\ where h is 
Planck's constant and v is the frequency of the photon, 
we have 
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We absorbed the h into the free parameters governing 
the diffusion and drift. We see that the rates of diffusion 
and drift in polarization angle are frequency dependent. 
Note that the Lorentz symmetry we have assumed is only 
invariance under the proper orthochronous component of 
the full Lorentz Group. The drift term explicitly breaks 
parity invariance. 

Our model has assumed that spacetime is Minkowski 
spacetime. Following [3] we can model the effect of an 
expanding universe by setting the frequency to depend on 
time, V = v{t), such that a{t)y{t) = a^vo where a{t) is 
the scale factor of the universe, oq is the current value of 
a and vq is the current (observed) value of the frequency 
of the photon. If we define a new time coordinate t' by 
dt'/dt = a{t)/aQ then our diffusion equation keeps the 
same form as ( 11 1 



dp 
dtf 



P- 



(12) 



For a matter dominated universe a ~ t2/3 and a range of 
t of 10^° Planck times becomes a range of t' of 3/5 x 10^° 
Planck times. We drop the subscript from v and the 
prime from t in what follows. 
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Consider a beam of photons of frequency v whose 
polarization is initially described by Stokes parameters 
(U,Q). The drift will result, after a time t in a beam 
whose polarization angle, $, has rotated by x •= td/v: 



$ = tan ^ 



= $ + X . 



(13) 



The diffusion will result in a decrease in the magnitude 
of polarization by a factor exp(— /i): 



(14) 



where fi can be calculated to be Atc/v. 

Clearly, such an effect would be most pronounced in 
photons that have propagated over long distances. Thus 
CMB photons, whose polarization is correlated over in- 
coming directions and have traveled over cosmological 
distances without rescattering, offer the best chance to 
constrain the parameters in the model. The polarization 
of the photons is imprinted during the last stages of re- 
combination as the photons decouple from baryons and 
enter the free streaming regime. The correlation in the 
polarization of photons arriving from different directions 
is encoded in a set of angular power spectra where 
XY are the different spectral cross-correlation in total in- 
tensity T and grad-type and curl-like components E and 
B respectively. The spectra are calculated by solving the 
full Einstein-Boltzmann system describing the evolution 
of perturbed fluids in a particular cosmological model 
[llj . The rotation and suppression of polarization along 
the trajectory will modify the angular power spectra of 
the observed CMB photons Cf^ — )■ ■ The mapping 
for each spectrum is given by [HI [13] 



(15) 
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= e 
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e-^^Cf^ sin(4x). 



where we have assumed no BB contribution to the origi- 
nal spectra (no primordial gravitational waves). The TT 
spectrum is not modified as it is not sensitive to the po- 
larization of the photons. 

A number of assumptions are implicit in the simple 
mapping given in ( [T5| ). Firstly it assumes that the pic- 
ture of recombination (when the polarization is imprinted 
on the CMB) and reionisation (a further source of po- 
larization) is not altered in this model. It assumes the 
background cosmological evolution is the same for a given 
set of cosmological parameters. In fact, polarization on 
large scales is also generated after the universe is reion- 
ized and this could introduce a mild scale dependence of 
the diffusion-rotation effect. 

To constrain this scenario with available CMB data 
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FIG. 1: The marginalized 2-d posterior density in the po- 
larization rotation angle xiso (in degrees) and polarization 
depth /ii5o for the reference frequency of 150 GHz. The con- 
tours indicate the area bounding 68% and 95% of the den- 
sity. We show the result for two data combinations; the first 
includes all polarization data and the second excludes the 
QUaD results. The results are consistent with no effect, how- 
ever including the frequency dependence reduces the tension 
between the QuaD and 'no QUaD' combinations highlighted 
in EOl. 



we modify the CosmoMC^ Monte Carlo Markov Chain 
(MCMC) package to fit for standard ACDM model pa- 
rameters together with polarization rotation x ^^nd po- 
larization depth fi. The standard parameters are: cold 
dark matter and baryonic matter physical densities flch^ 
and flbh^, angular diameter distance measure 0, optical 
depth to reionisation r, and primordial scalar perturba- 
tion amplitude As and spectral index Tig. We assume a 
uniform prior of sufficient range in each parameter and 
do not include any primordial tensor contributions. How- 
ever we do fit to TB, EB, and BB data since these are 
not expected to vanish any longer in the modified model. 

We fit to a combination of data which includes all po- 
larization sensitive experiments which have reported a 
detection of the EE power. These are the DASI results 
[Tl] , the final CBIpol rcsuhs [H], the Boomerang 2003 
flight results [16, , the WMAP 5-year results [H, and the 
latest BICEP [H] and QUaD [11] results. Of these, the 
last two contain the highest signal-to-noise determina- 
tion of the polarization spectra on scales below a degree. 
Both BICEP and QUaD have published TB and EB data 
which are crucial in constraining polarization rotation ef- 
fects [I8l[20] . We include the published TB and EB band 
powers with band power window functions mimicking the 
published TE and EE ones. The frequency dependence 
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TABLE I: Marginalized 1-d constraints for the polarization 
rotation angle x ^'^'^ 95% upper limits for the polarization 
depth II. 





ALL Pol 


No QuaD 


Pol rotation: x (degrees) 
Pol depth: fj, 


-o.o5t°:f3 

< 0.024(95%) 


1 11+0.55 
~-L--LJ^-0.55 

< 0.026(95%) 



is accounted for by scaling the effect for each experiment 
to a reference frequency of 150 GHz. 

The MCMC chains sample the posterior density in the 
8-dimensional parameter space. Once the sampling has 
converged we integrate the densities over the standard 
parameters which gives the marginalized posterior in xiso 
and /ii5o. We show the result of the marginalization in 
the Xi50 Mi50 plane in Fig. [T] We plot the 68% and 
95% density contours for two data combinations. The 
tightest constraints are obtained by the combination of 
all polarized data and are consistent with no rotation and 
vanishing polarization depth. We also include the results 
for the case where the QUaD data is excluded. Although 
the 'no QUaD' result, driven mainly by the BICEP and 
WMAP measurements, prefers a non-zero rotation an- 
gle, the indication is weaker than that reported in |20j . 
However we do recover their result when not accounting 
for the frequency dependence. This may be an indication 
of a frequency dependent effect in the data whereby the 
lower frequency WMAP data tend favour less rotation 
given its effect is roughly twice as large compared to the 
reference frequency of 150GHz. Future observations at 
multiple frequencies will easily determine whether this is 
the case. 

Table[T]shows the 1-d marginalized constraints on Xi50 
and /xi5o showing the marginal bias towards a negative 
rotation driven mainly by the BICEP and WMAP data. 
The "All polarized" result is consistent with no effect. 

The model building strategy employed here is based 
on the assumption that the fluctuations in dynamical 
variables due to spacetime uncertainty are small enough 
that they result, in the hydrodynamic approximation, in 
a continuous, Brownian motion through the state space. 
If spacetime uncertainty causes more violent, discontinu- 
ous jumps in the variables, this will have to be modeled 
by a Boltzmann equation rather than a diffusion equa- 
tion. For this reason and also as data begins to constrain 
our phenomenology, the need for microscopic models that 
can give us a handle on the parameters from more fun- 
damental physics becomes more acute. For example, an 
improvement on point particle models would be to treat 
particles as wave packets of a scalar field on a causal set 
using the discrete D'Alembertian operators described in 
[H1[71[5T]. Assuming continuity, however, the power of 
our model is its robustness: if photons can be modeled 
as classical particles with polarization degrees of freedom, 
then any Lorentz invariant effect of underlying spacetime 
uncertainty - whether due to discreteness, fluctuations. 



fuzziness, foaminess or whatever - will be of this form. 
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